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Anderson's orthogonality catastrophe in graphene, at energies close to the Dirac point, is analyzed. 
It is shown that, in clean systems, the orthogonality catastrophe is suppressed, due to the vanishing 
density of states at the Dirac point. In the presence of preexisting localized states at the Dirac 
energy, the orthogonality catastrophe shows similar features to those found in normal metals with a 
finite density of states at the Fermi level. The implications for the Kondo effect induced by magnetic 
impurities, and for the Fermi edge singularities in tunneling processes are also discussed. 
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INTRODUCTION. 

Graphene has attracted a great deal of attention re- 
cently, due to its novel fundamental properties and po- 
tential applications P, [3, H, |3j 0|- It is by now well es- 
tablished that its electronic properties at low energies 
are well described by the two dimensional Dirac equa- 
tion. At half filling graphene should be a semimetal, 
with a vanishing density of states. This fact implies that 
many properties of a metal, which are parametrized by 
the density of states at the Fermi level, are different in 
a clean graphene sample. The description of the elec- 
tronic bands in graphene based on the Dirac equation 
also leads to localized states in samples with edges 0] 
or lattice defects [1, 0. These states change the density 
of states near the Dirac energy, as they induce a peak at 
this energy. Hence, the density of states of graphene at 
the Dirac energy can either vanish, in a clean sample, or 
diverge, if localized states are induced.. 

We study here Anderson's orthogonality catastrophe 
fAOC)[loj in clean and dirty graphene. The AOC can 
be considered the simplest non trivial feature in the re- 
sponse of a metal, and it is dependent on the value of 
the density of states at the Fermi level. The AOC di- 
rectly leads to many singularities in experiments which 
probe the dynamical response of a metal, like the Fermi 
edge singularity in X-ray absorption [ill and singu- 
larities in the transport properties of quantum dots and 
metallic grains[l3, 13, EE]- In graphene, the interplay 
between the AOC and Coulomb blockade may be rele- 
vant for the analysis of transport experiments on small 
quantum dots 0, [3l ■ 

The Kondo effect induced by magnetic impurities in 
metals can be seen as a direct consequence of Anderson's 
orthogonality catastrophe [l7|, ll8|. The coupling between 
the impurity spin and the conduction electrons can be di- 
vided into a transverse term, J_l, which leads to spin- flip 
processes, and a longitudinal term, J|| , which induces an 
AOC associated to the same spin flips. This AOC leads 



to a strong suppression of spin fluctuations, although the 
effects of J± prevail at the lowest temperatures. These 
two competing processes can be defined, in a very trans- 
parent way, in the dissipative two level system [l9|, which 
is equivalent to the Kondo Hamiltonian. The Kondo tem- 
perature, Tk, can be seen as the scale at which spin flip 
processes ultimately cut off the AOC. 

The AOC is modified in disordered metals [20l| and 
ballistic mesoscopic systems, due to the changes in the 
electronic wavefunctions 21 , . We will analyze the 



AOC i n g raphene u sing the numerical methods explained 
in 2]J, 122] (see alsoj23| ) , and also a phaseshift analysis 



similar to that in|lC 

We analyze first the phaseshifts induced by a local po- 
tential, first in clean graphene, and then in graphene in 
the presence of preexisting localized levels. The next sec- 
tion presents a numerical study of the full overlap be- 
tween the electronic ground state before and after the 
potential is turned on, and an analysis of the scaling of 
this quantity with system size. The last section discusses 
the main implications of our work for the Kondo effect 
in graphene, and the Fermi edge singularities associated 
to tunneling processes. 



PHASESHIFT ANALYSIS. 

The overlap S between the Slater determinants which 
describe the electronic wavefunction before and after a 
potential is switched on can be written as[lo|: 



(1) 



where N is the number of electrons, and 5i is the phase- 
shift induced by the potential in the scattered waves at 
the Fermi level with angular momentum I. In a typical 
metal, a weak local potential of strength eg ^ £f induces 
a phaseshift in the s channel which can be approximated 
by 6o « eoN{eF) ^ 1, where N{eF) is the density of 
states at the Fermi level. 
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This analysis can be extended in a straightforward way 
to graphene, where the electronic wavefunctions can be 
approximated by the two dimensional Dirac equation (see 
below): 

where the two signs correspond to the two inequivalent 
corners of the Brillouin Zone of the honeycomb lattice. 

We use eq.([T]) in order to describe the dependence of 
the overlap on the number of electrons, by computing 
analytically the phaseshifts induced by different types of 
potentials. In the following, we use energy and momen- 
tum units such that wp = 1. 



Clean graphene. 

We analyze first the phaseshifts induced by a circular 
potential well in clean graphene, and we describe the elec- 
tronic wavefunctions using the continuum Dirac equa- 
tion, eq.([2) . We expect that this approximation will 
describe qualitatively the effects of a local perturbation 
in the graphene lattice. 

We assume that the potential well can scatter elec- 
trons between the K and K' valleys, as it is the case for 
sufhciently localized potentials in graphene. 

Using cylindrical coordinates, the Hamiltonian in the 
clean system can be written as: 
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where the two first entries correspond to the K point, 
and the two last ones to the K' point. 

We add a constant perturbation in the region r < Rq: 
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where eo is a constant energy shift, and A is a potential 
which induces scattering between the two valleys, and 
it is c omp atible with the symmetries of the honeycomb 
lattice^^. 

We analyze the scattering of an incident s wave with 
incoming energy k: 



*inc(r, (f>) 
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(5) 



where Jo{x) and Ji{x) are Bessel functions of the 



first kind. They satisfy: lim^j^o Jo{x) 
linij^^o Ji{x) ~ x/2. 



1, and 



The reflected waves outside the well can be written as: 
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Yo{x) and Yi{x) are Bessel functions of the second kind. 
They satisfy: lim^j-^o ^o(a^) ~ 2/7r(log(x/2) -f 7), and 
\mix—>QYi(x) « —2/(Trx). The first contribution on the 
right hand side of eq. ([H]) is a reflected wave in the same 
valley, and the second term is a wave in the opposite 
valley as the incident wave. 

Inside the potential well, the spectrum has a gap for 
energies eo — A < e < eo-l-A. Within this range of 
energies, the wavefunction inside the well can be written 
as: 
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Io{x) and7i(a;) are modified Bessel functions of the first kind. They satisfy: lim^j^o ^o(a^) ~ Ij and limj;^o -^i (2;) ~ x/2. 
The value of k' in eq. ([7]) is given by: e = VA^ — fc'^. As fc = e + eo , we have k' — -y/A^ — (fc — eo)^- 
For |e — eol > A, we have: 



*trans(»', ' 
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(8) 



and e = VA2 + fc'^, and fc' = ^(fc - eo)2 - A2. 

The scattering phaseshifts are determined by the reflection coefficients Ri and R2 defined in eq. © . The boundary 
conditions at r — Rq are simply the continuity of the spinors, which define a set of four equations for the four variables 
i?i,i?2,Ti and T2. 



For A = we have i?2 = ^2 = and i?i — R. As lini2;^oo Jo{x) ~ ■y/2/(7rx) cos(a; — 7r/4), and lima;^oo i1](a;) 
;he phaseshift (5 is t 

tan((5) = i?(fci?o) = 



y^2/{TTx) sin(x — 7r/4), the phaseshift (5 is tan((5) = R. We find 

Ji[(fc - eo)Rn]MkRo) - Jo[{k - en)Ro]Ji{kRo 



Ji[{k ~ eo)Ro]Yo{kRo) ~ Jo[(fc - eo)i?o]n(fci?o) 



(9) 




Energy (kRJ 

FIG. 1: (Color online). Reflection coeflicients of a circular 
well with eoRo = 0.5 and ARo = 0.1. 



and: 

lim R{kRo) = tan(eoi?o) (10) 

kRo — ^00 

Results for eoRo ~ 0.5 and ARq = 0.1 are shown in 
Fig.p. In all cases, with or without (eq.®) interval- 
ley scattering, the reflection coefficients vanish at the 
Dirac point, fc = 0. This result can be simply under- 
stood by noting that a finite refiection coefficient implies 
a reflected wavefuntion with a component Yi(kr) which 
diverges as fc 0. The phaseshift vanishes linearly as 
fc — > 0, in agreement with general arguments based on 
the vanishing of the density of states at the Dirac point. 



The vanishing of the phaseshift at the Dirac point im- 
plies that the overlap between the Slater determinants 
before and after the potential is switched on does not 
scale like some power of the number electrons, and the 
AOC does not take place at this energy. 

Phaseshift analysis in the presence of a locahzed 
state. 

We will neglect here possible intervalley scattering 
terms. We study the phaseshifts induced by a weak po- 
tential near the edges of a circular void which supports 
surface states. A sketch of the model is shown in Fig. [2]. 
We write the wavefunction as: 

The edge of a crack, or extended vacancy is modeled by 
the boundary condition: 

i/'i(r) = 0, ren (12) 

where is the boundary of the void. 

We analyze a circular void, of radius R' . The boundary 
condition, eq. (fT2|) . allows for solutions at zero energy of 
the type: 




where the two signs correspond to the two inequivalcnt 
corners of the Brillouin zone. 

Eq. fT^ implies, for s-wave scattering: 

aoJo{kR')+l3QYoikR')=0 (14) 
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The phaseshift induced by the void, before the potential 
whose effect we want to calculate is turned on, is: 
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FIG. 2; (Color online). Sketch of the model with circular 
symmetry used to study the AOC in the presence of localized 
levels. An infinite potential exists for < r < R' , mimicking a 
vacancy. The perturbation leading to the AOC is modeled as 
a constant potential, eo, for R' < r < R (see text for details). 



Next, we model a weak impurity near the void as an 
isotropic perturbation of depth eo, defined in the region 
R' < |f I < R- Following eqs. ()5l6p and neglecting inter- 
valley scattering, the wavefunction can be written as: 



a'M{k + eo)r]+(3'Yo[{k + eo)r] \ . . o 

a'Ji[(fc + eo)r]e'* + /3'yi[(fc + eo)r]e'V " " . . 

aJoikr) + PYoikr) \ ^'"'^ 

aJi (fcr)e''^ + /3Yi (fcr-)e^'^ ) - ' 

with boundary conditions: 

a'Jo[(fc + eo)i?']+/3'Fo[(fc + eo)i?'] = 
a'Jo[{k + to)R] + P'Yo[{k + eo)R] = aJo{kR) + pYo{kR) 

a'Ji[(k + tQ)R]+l3'Yi[{k + €o)R] = aJi{kR) + f3Yi{kR) (17) 




These equations allow us to obtain the phaseshift of 
the combined system, void and circular impurity, as 
S = arctan(/3/Q;) . The overlap between the Slater de- 
terminants before and after the impurity potential is 
switched on, is determined by the phase difference, S—Sq, 
where 6a is given in eq. ljlSp . 

Results for the individual phaseshifts 6 and 6q, as well 
as their difference are shown in Fig. [5] for eo = 0.1, i?' = 
0.9 and R = 1. In this regime of energies much lower 
than eo, the phaseshift 5 seems to approach i5o from be- 
low, indicating that the repulsive character of the void is 
weakened by the additional constant potential. For the 
small energies close to the Dirac point focused on here, 
the relative phaseshift, 6 — Sq, is always finite and seems 
to approach a constant. This behavior differs strikingly 
from our findings for clean graphene where the vanishing 
of the phaseshift at the Dirac point (cf. Fig.[T]) indi- 
cates the suppression of AOC. In the presence of voids, 
the small dependence of the phaseshift induced by an 
additional external potential on energy near the Dirac 
point implies that the overlap between Slater determi- 
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FIG. 3: (Color online). Phaseshift (5o induced by a void, the 
phaseshift delta resulting from the additional switching on 
of a constant potential, and the resulting relative phaseshift 
(5 — (5o induced by a circular impurity potential surrounding a 
void (see text for details). 
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Overlap (12 x 12 cluster) 




FIG. 4: (Color online). Overlap and effective phaseshift as 
function of filling and potential strength (see text for details). 



nants should scale with the number of electrons in a 
similar fashion to that in a normal metal with a finite 
density of states. We shall see in the remainder of this 
paper that there are indeed considerable differences be- 
tween clean graphene vs. graphene with localized states, 
that are visible, e.g., in the behavior of the AOC overlap. 



CALCULATION OF THE OVERLAP. 
Clean graphene. 

The overlap between the unperturbed and perturbed 
Slater determinants for clean graphene clusters of differ- 
ent sizes has been calculated using the methods described 
in 21, 22, 2^. The perturbation is a local potential at a 
given site, A = cq. Its strength is measured in terms of 



the scaled perturbation strength (x A/d with d being the 
mean level spacing 6/[(N (N-|-l)-2]. Periodic boundary 
conditions are used in systems with N x N unit cells, 
up to = 80; the vertical stripes visible in Fig. [4] are 
an artefact of the periodic boundary conditions. The re- 
sults for the overlap for iV = 12 and different potential 
strengths (ranging from weak to strong for repulsive as 
well as attractive perturbations) are shown in Fig. [4]. An 
effective phaseshift can be defined by dividing the energy 
shift of the level closest to the Fermi energy by the av- 
erage level spacing in that energy range. This phaseshift 
is also shown in Fig. [3]. 



FIG. 5: (Color online). Scaling of the overlap with cluster 
size N at the Dirac point (DP, ep — 0, corresponding to a 
filling of 0.5, top panel) and away from the Dirac point (at 
fixed filling 0.55, corresponding to ep ~ 0.67, lower panel). 
See text for details. 



The dependence of the overlap with system size is dif- 
ferent at the Dirac point from that at other energies. 
This dependence is shown in Fig. [5]. The overlap is 
almost independent of system size at the Dirac point, 
cf. the upper panel. This result is consistent with the 
phaseshift analysis, which shows that the phaseshift van- 
ishes at the Dirac point. Indeed, AOC is suppressed at 
the Dirac point. Away from the Dirac point, the con- 
ventional behavior of the AOC overlap is recovered, see 
the lower panel of Fig. [5]. To this end, AOC overlaps 
for fillings ranging from 0.54 to 0.56 were averaged over. 
Clearly, the AOC overlap is no longer suppressed and ap- 
proaches zero in the thermodynamic limit following the 
well-known power-law dependence on the number of par- 
ticles (cx [N{N + 1)- 2]), cf. inset of Fig.[S]. 



Graphene with localized states. 

The method described in [2^ assumes that the wave- 
functions of all eigenstates of the unperturbed system 
have the same weight on the site where the perturbation 
is turned on. This leads to a considerable simplification 
of the calculation of the overlap between Slater deter- 
minants. Generalization of this method generalized to 
to chaotic mesoscopic systems 2lJ, |22| was done based on 
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FIG. 6: (Color online). Dependence of the overlap on pertur- 
bation strength when the perturbation is turned on near an 
existing vacancy (empty circles, black), and in clean graphene 
(filled circles, red) . Calculations are done for 12 x 12 clusters. 
Circles correspond to one hole in the cluster (Dirac energy, 
ep = 0), whereas diamonds characterize a cluster with five 
holes (corresponding to ep = —0.5, or a filling of ~ 0.47). 



FIG. 7: (Color online). Dependence of the overlap on pertur- 
bation strength, at the Dirac energy, when the perturbation 
is turned on near an existing vacancy (circles, black), and in 
clean graphene (diamonds, red). Large symbols correspond 
to a 12 X 12 cluster and small symbols correspond to a 15 x 15 
cluster. 



CONCLUSIONS. 



the statistical properties of the chaotic wave functions. 

In the presence of a defect which induces a localized 
state, like a vacancy, the wavefunctions of the unper- 
turbed state, where the localized state is already present, 
do not possess translational symmetry. Therefore direct 
diagonalization and calculation of the overlap determi- 
nants was used for the study of clusters of moderate sizes. 



Results for the overlap for clusters with 12 x 12 unit 
cells are shown in Fig. [6]. At the Dirac point, the pres- 
ence of a vacancy, which induces a localized state, en- 
hances significantly the dependence of the overlap on the 
strength of the potential. Away from the Dirac point, 
the difference in the overlap with and without a vacancy 
is much less pronounced. We show in Fig. [7] the depen- 
dence of the overlap with cluster size, at the Dirac energy. 
As anticipated in the discussion of Fig. [3] , the presence of 
a vacancy near the potential which is turned on modifies 
significantly the results in comparison with a clean sys- 
tem. In the latter, the dependence on size is negligible, 
in agreement with the results shown in Fig. [5]. There is, 
on the other hand, a substantial dependence on cluster 
size when a vacancy induces a localized state at the Dirac 
energy. 



The results presented here show the existence of 
two regimes for Anderson's orthogonality catastrophe in 
graphene at low fillings, depending on whether there are 
localized states at the Dirac energy or not. In the ab- 
sence of localized states the AOC is suppressed near the 
Dirac point, in agreement with the vanishing of the den- 
sity of states at this energy. When localized states are 
present, the AOC is qualitatively similar to that found in 
metals with a finite density of states. The latter behavior 
is a consequence of the fact that, when localized states 
are sufficiently near the Fermi surface, they contribute 
to the non adiabatic response of the electron gas. This 
situation is unique to graphene, as, in most metallic sys- 
tems, localized states appear at energies well below the 
Fermi level. 

The features discussed above imply that the Kondo 
effect in graphene also depends on the strength of the 
scalar potential induced by the magnetic impurity. If 
the potential induced on the graphene electrons is weak, 
as when the magnetic impurity is at some distance of 
the graphene plane, we expect the formation of a Kondo 
resonance to be suppressed, and the magnetic impurity 
will give rise to a free magnetic moment. On the other 
hand, if the magnetic impurity lies within the graphene 
plane, it will give rise to a strong scalar potential, and 
possibly to localized states at the Dirac energy. Then, 
the Kondo effect will not be suppressed, despite the low 
density of states in graphene near the Dirac energy. 
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Similar effects can be expected for the Fermi edge sin- 
gularities induced by electrons tunneling into or out of 
graphene quantum dots. The strength of the Fermi edge 
singularities depend on the existence of localized states 
in the quantum dot. These states will be induced in 
graphene dots with sharp and rough edges, where, in 
addition to Coulomb blockade, the AOC associated to 
electron tunneling will further suppress the conductance 



at low voltages [13|,[1 
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